Exponential Levy-type models with stochastic volatility and 



stochastic jump intensity 



Matthew Lorig * 
This version: March 3, 2013 



Abstract 

We consider the problem of valuing a European option written on an asset whose dynamics are de- 
scribed by an exponential Levy-type model. Both the volatility and jump-intensity are allowed to vary 
stochastically in time through common driving factors - one fa st-varying and one slow-varying. Our re- 



sults e xtend the class of multiscale stochastic volatility models of 



(2011 



Fouaue. Papanicolaou. Sircar, and Solna 



) to models of the exponential Levy type. Using generalized Fourier transform techniques we de- 
rive an explicit formula for the approximate price of a European-style derivative. For derivatives with 
smooth and bounded payoffs we establish the accuracy of our pricing approximation. As an example 



Mertonl (|l976t ) to include stochastic volatility 



of our framework, we extend the jump-diffusion model of [ 
and stochastic jump-intensity. We perform a calibration to S&P500 options and find that the extended 
Merton model with a single fast-varying factor of volatility provides a better fit to implied volatility 
than both the traditional Merton model and the fast mean-reverting stochastic volatility models of 

( 201lh . 



Fouaue. Papanicolaou. Sircar, and Solna 
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1 Introduction 



An exponential Levy model is an equity model in which an underlying S — e x is described by the exponential 



of a L evy process X. Such models extend the geometric Brownian motion description of 



Black and Scholes 



(see, 



([1973 ) by allowing the underlying S to experience jumps, the need for which is well-documented in literature 



Erakerl (|2004f ) and references therein). In addition to allowing the underlying S to jump, exponential 
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Levy models are important because they capture many of the stylized features of asset prices, such as heavy 
tails, high-kurtosis and asymmetry of log returns. 

Several well-known models fit within the expo nential Levy class: the jump- diffusion model of 



1976), the pure jump models of 



Mandelbrot J 1963 ), the variance gamma model o f 



19981 ). the e xtended Ko p onen family o f 



tial model of 



Koul (120021) . 



Lewis! ( 2001 ) 



Morton 



Madan. Carr. and Chang 



Bovarchenko 



Lipton 



and Levendorskiil (|2000l ) and the double exponen- 



20021 ) show that all of the above-mentioned models al- 



low for fast and easy computation of European option prices via one-dimensional Fourier transforms. A 
comprehensive reference on the subj ect of option-pricing in a n exponential Levy setti ng can be found in 



Bovarchenko and Levendorskiil (|2002l ). as well as Chapter 11 of ICont and Tankovl (|2004l ). 



Despite their success, exponential Levy models have some shortcomings. For example, because the 
log returns of any exponential Levy process are independent and identically distributed, these models 
cannot exhibit volatility clustering (the tendency for volatility to rise sharply for short periods of time) 
or the leverage effect (the tendency for volatility to rise as asset prices fall); both of these phenomena 
are well-documented in literature. Levy processes also exhibit constant jump intensities. However, a 
recent study of S&P500 index returns i ndica tes that jump-intensities, like vol atility, are sto c hastic (see 



Christoffersen. Jacobs, and Ornthanalail (j2009f ) ) . To address these shortcomings. ICarr and Wul (|2004l ) add 



stochastic volatility (with corr elation to the underly ing) by stochastically time-changing a Levy process. No- 



tably, the models described in ICarr and Wul (|2004f) maintain the analytic tractability that makes the class 
of exponential Levy processes attractive. 

In this paper, we address the need for volatility clustering, the leverage effect and stochastic jump 
intensity by modeling the returns process X by a Levy- type process whose local characteristics (7t, <xt, z^t) 
are stochastic. We then use generalized Fourier transform techniques and singular perturbation methods to 
derive an explicit formula for the approximate price of a European-style derivative. 

Much like geometric Brownian motion arises as special case of an expon ential Levy process, the class of fast 



mean- rever ting and multiscale stochastic volatility models considered in 



pOOOh and 



Fouaue. Papanicolaou, and Sircar 



Fouaue. Papanicolaou. Sircar, and Solnal (|201ll) arise as a special subset of the class of mod- 



els we consider. In fact, by removing jumps from our fr amework, one recovers the Fouri e r repr esen- 



tation of the European option pricing formul as derived in 



Fouque. Papanicolaou. Sircar, and Solna 



Fouaue. Papanicolaou, and Sircarl (|2000l ) and 



(|201ll) . 

The rest of this paper proceeds as follows. In section [2] we introduce a class of exponential Levy-type 
models in which the volatility and jump-intensity are stochastically driven by a common fast- varying factor. 
In section [3] we derive an expression for the approximate price of a European option (Theorem 13- L[) when 
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the underlying is described by the class of models introduced in section [2] We also quantify the accuracy 
of our pricing approxim ation (Theorem 13.21) . In section 2J as an example of our framework, we extend the 



jump-diffusion model of 



Mertonl (|l976l ) to include stochastic volatility and stochastic jump-intensity. We 



also compute (numerically) the implied volatility surface generated by this example. A calibration of the 
extended Merton model to S&P500 index options is performed in section [5) In section [6] we briefly describe 
how the class of models introduced in section [5] can be extended to allow for multiple driving factors of 
volatility and jump-intensity - one fast-varying factor and one slow-varying factor. Proofs are provided in 
an appendix. 



2 Stochastic volatility and jump intensity Levy-type processes 

Let (f2,£F, P) be a probability space endowed with a filtration F = {£ft,f > 0}, which satisfies the usual 
conditions. Here, P is the risk-neutral pricing measure, which we assume is chosen by the market. The 
filtration F represents the history of the market. For simplicity, we assume that the risk-free rate of interest 
is zero so that all non-dividend paying assets are (P, F)-martingales. All of our results can easily be extended 
to include constant or deterministic interest rates. 

We consider a non-dividend paying asset S whose dynamics under P are described by the following 
Ito-Levy stochastic differential equation (SDE) 

dS t = a(Yt)S t dW t + S t - I (e z - 1) dN t {Y t ,dz), S = logx, 

(under P) 



\a{Y t ) - -A(Y" t ) P(Y t ) ) dt + -f3(Y t )dB t , Y Q = y, 



(2-1) dY t -- 

d(W,B)=pdt, |p|<l- 
Here W and B are correlated Brownian motions and N(Y, dz) is a compensated Poisson random measure 

dN t (Y u dz) = dN t (Y t ,dz) ~ C(Y t )v(dz)dt, E[dN t (Y t ,dz)\Y t ] = ((Y t )is(dz)dt, 

We require that the measure v satisfy 

min(l, z 2 )v(dz) < oo, / e z v(dz)<oo, and / \z\v{dz) < co. 

I J\z\>l ^|z|>l 

The first integrability condition must be satisfied by all Levy measures. The second integrability condition is 
needed to ensure E[St) < oo for all t £ R + . The last integrability condition allows us to replace the indicator 
function that usually appears in the Levy-Kintchine formula I{| z |<i} with the constant 1. Although we do 
not require it, a correlation of p < between W and B would be consistent with the leverage effect (i.e. a 
drop in the value of S will usually be accompanied by an increase in volatility). 
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Note that both the volatility of S, given by <j(Y), and the jump-intensity ((Y)i/(dz), which controls the 
jumps of S, are driven by a common stochastic process Y. We require the existence of a constants a, a and 
C such that < a < a(y) <W < oo and < ((y) < ( < oo for all y in the state space of Y. 

2.1 The driving process Y and model assumptions 

The process Y is fast-varying in the following sense. Under the physical measure P, the dynamics of Y are 
described by 

dY t = \a(Y t )dt+-l3(Y t )dB t } (under P) 
where B t = B t — J Q A(Y s )ds is a P-Brownian motion. The generator of Y under P is scaled by a factor of 

A{. = ±(±l3 2 (y)d 2 yy + a( y )d y y 

Thus, Y operates with an intrinsic time-scale e 2 . We assume e 2 << 1 so that the intrinsic time-scale of Y is 
small. Thus, Y is fast-varying. Let Y 1 be a diffusion process whose infinitesimal generator is Ay (so that, 
in distribution, Y t = Y^, 2 under P). We assume the following: 

1. Under P, the process Y 1 is ergodic and has a unique invariant distribution Fy- 

2. The smallest non-zero eigenvalue of — A Y is strictly positive. 

3. There exists a constant C(k) < oo such that 

supE^Y] < C(k). 
t 

4. The coefficients a{y) and P(y) are smooth and at most polynomially growing in y. 
Typical processes which satisfy assumptions [U [2] and [3] above are 

OU process : a(y) = m — y, f3(y) — v, Fy ~ Normal, 

CIR process : a(y) — m — y, (5{y) = v^/y, Fy ~ Gamma. 

We further assume: 

5. The functions <r(y) and £(?/) are smooth and that the solutions rj(y) and t;(y) to Poisson equations 
(2.2) AoV = v 2 -(<r 2 ), A>e = C-(C), (f) ■= J f(v)Fy(dy). 

6. The market price of volatility risk A is smooth and bounded above and below. 
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7. The existence of a unique strong solution to SDE (|2.ip . 
As mentioned in the introduction, the class of models described by (|2.ip are a natural extension of those 



considered in 



we consider and those considered in 



l-'onquo. Piipauicolnon. and Sircai 1 200(1 . T he key dilFevences l>ot\v<vn |h<> class of models 



Fouaue. Papanicolaou, and Sircarl (j2000h are that (i) we allow for the 



underlying S to jump and (ii) we allow for the jump intensity to be stochastic. 

The use of a fast-varying driving factor Y is warranted when the life of an option t is large compared 
to the characteristic time-scale e 2 of the d r iving p rocessfi.e., t ^> e 2 ). A variogram analysis performed 



by iFouque. Papanicolaou. Sircar, and S0mal (|2003aJ) indicates a fast- varying component of volatility with a 
characteristic time-scale of 1.7 days. Typically, then, the use of models with fast-varying driving processes 
is limited to options with expiricies of greater than 17 days. 



3 Option pricing 

We wish to price a European-style option, which pays H(St) at the maturity date t > 0. It will be 
convenient to introduce th e returns process X = log S. Using Ito's formula for Ito-Levy processes (see 
0ksendal and Suleml (|2005l ). Theorem 1.14) one derives 



dX t = j(Y t ) dt + a{Y t ) dW t + z dN t (Y t ,dz) 



X = x, 



where the drift j(Y t ) is given by 



l(Yt) = — e\Yt) C(Yt) (e z -l- z)v{dz). 
Using risk- neutral pricing, the value u e {t, x, y) of the European option under consideration is 



l {t,x,y) =E X>2/ [h{X t )\ 



h{x) :=H{e x ). 



From the Kolmogorov backward equation we find that u e (t,x,y) satisfies the following partial integro- 
differential equation (PIDE) and boundary condition (BC) 



(3.1) 



(-d t +A e )u E (t,x,y) = 0, 



'(0,x,y) = h(x). 



Here, the partial integro-differential operator A e is the generator of (X,Y), which is defined on functions 
/ : M 2 -> E by 



A s f(x,y) := lim \ (E XtV [f(X u Y t )] - f(x,yfj 



(if the limit exists) . 
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If / € Cq(M 2 ), the space o f twice differentiable f unctio ns with compact support, then the above limit does 
exist (see Theorem 1.22 of 0k sendal and Suleml ((2005)) and A e is given by 



(3.2) A e = -^A + -A! + A 2 , 

Ao=A 1 Y = ^ 2 (y)d 2 yy + a(y)d v , 
A 1 =p/3(y)a(y)d 2 xy -A(y)(3(y)d y , 

(3.3) A 2 = 1 {y)d x + \<7 2 {y)d 2 xx + ((y) [ (6,-1- zd x ) u(dz). 

z JR 

Here, 9 Z is the shift operator, which acts on functions / G dom(A E ) via 9 z f(x,y) := f(x + z,y). From now 
on, we define A £ f(x,y) by expressions (|3.2[) - (13.3[) . for all / such that the partial derivatives and integrals 
exist at (x, y). 

To establish the accuracy of our pricing approximation (Theorem l3.2|) we need for Cauchy problem (|3.ip 
to have a unique point-wise classical solution. As such, we shall assume enough regularity of the coefficients 
(a, (3, A, a, £) and the payoff function h so as to ensure that (|3.1[) has classical solution. This will be the case, 
for example, if the coefficients /3 2 , a, /3a, A/3, 7, a 2 , £, the payoff function h and its first two derivatives 
are bounded. Note that the above conditions are sufficient, but not necessary, to ensure a classical solution. 
Please see Theorem IB. II of Appendix [B] 

3.1 Formal asymptotic analysis 

For general (a, (, a, (3, A) there is no analytic solution to (|3.ip . We notice, however, that terms containing 
e in (|3.1I) are diverging in the small-e limit, giving rise to a singular perturbation about the 0(1) operator 
-dt +^.2)- This special form suggests t hat we seek an asymptotic solution to PIDE (|3.1[) . Thus, following 



Fouque. Papanicolaou, and Sircar 



20001 ) . we expand u e 



(3.4) u E = 



u r , 

71=0 

Our goal will be to find an approximation u e = uq + e u\ + 0(e 2 ). The choice of expanding in integer powers 
of e is natural given the form of A E . We will justify this expansion when we prove the accuracy of our pricing 
approximation in Theorem [321 

In the formal asymptotic analysis that follows, we insert expansion (|3.4[) into PIDE (|3.1|l and collect 
terms of like powers of e, starting at the lowest order. The 0(l/e 2 ) and 0(l/e) terms are 

0(l/e 2 ): Q=A n u , 
0(l/e) : 0=A 1 u o +A u 1 . 



G 



Noting that all terms in Aq and A\ take derivatives with respect to y, we choose uq = uo(t,x) and u\ = 
U\(t,x). Continuing the asymptotic analysis, the 0(1) and 0(e) terms are 

(3.5) 0(1) : = {-d t + A 2 )u + A a u 2 , 

(3.6) 0(e) : = (-d t + A 2 ) Ul + A lU2 + A u 3 , 

where we have used the fact that Ami = in the 0(1) equation. Equations (|3.5|) and (|3.6|) are equations of 
the form 

(3.7) A Q u = x- 

A solution u to (|3.7I) exists if and only if \ satisfies the centering condition 

(x) ■= J xdF Y = 0. 

For th e origin of the centering condition we refer the reader to equations (3.12)-(3.13) in 



Fouaue. Papanicolaou. Sircar, and Sol 



(|201ll ) and the discussion therein. Applying the centering condition to (|3.5[) and (|3.6|) yields 



(3.8) 0(1): = (-d t + (A 2 ))u , 

(3.9) 0(e): = (-8 t + {A 2 )) Ul + {A x u 2 ). 



Note, from and (|3.5|) and (|3 . 8|) we have 



A u 2 = ~(-d t +A 2 )u + (-d t + (A 2 ))u = -(A 2 - (A 2 ))u 
= -\(o- 2 -(° 2 )){dl x -d x )u. 



- (C - (C» (- J r (f - 1 - *) u(dz)d x +J(e x -1- zd x ) v(dz)^j u 

(3.10) = -Aq Qjj (d 2 xx - 8 X ) - $j (e z - 1 - z) u(dz)d x +tJjo,-l- zd x ) v(dzYj u 
where we have used r)(y) and £(y) as solutions to (|2.2|) . Thus, from (|3.9p and (|3.10p we find 

(3.11) 0(e) : (-dt + (A 2 )) Ul = -Su , 
where the operator 23, which is defined to act on functions of x (and not y), is given by 

23 = l-A x (d 2 xx -d x )-if (e* - 1 - z)u(dz)d x + £ / (e z - 1 - z9,)^(dz) 

(3.12) = V 3 (d 3 xxx - d 2 xx ) +U 3 (-J(e'-1- z)u(dz)d 2 xx + J (o z - 1 - ^^(dz)) 

+ y 2 (9L - d x ) +U 2 (- [ (e z - 1 - z) ^(cb)^ + / (t 



- 1 - Z&j ) l/(dz) , 



and the constants (V3, U3, V2, U 2 ) are defined as 



V a = -£(p*a v i 1 ), U 3 = -p((3ad y (i), V 2 = ±(pAd y r)), U i = (pAd v $. 



Note that equality (|3.12[) holds only when 25 is applied to functions that are independent of y. This is as far 
as we will take the asymptotic analysis. To review, we have found that Uo(t, x) and u\(t,x) satisfy PIDEs 
(|3.8j) and (|3.11j) respectively. We also impose the following BCs 



(3.13) 
(3.14) 



0(1) 

0(e) 



u o (0,x) = h(x), 
ui(0,x) = 0. 



3.2 Explicit solution for uo(t, x) and ui(t,x) 

In order to find explicit formulas for uo(t,x) and Ui(t, x), we note that the operator 

1 



(3.15) 



(^ 2 > = (7)& + §(a 2 )9L + (C> 



1 — zdx) v(dz). 



is the generator of a Levy process with Levy triplet ((7), (c 2 ), (C) v )- Thus, we may apply standard results 
from the classical theory of generalized Fourier transforms to obtain solutions to PIDEs (I3.8[) and (|3.11l) . 
Recall that the generalized Fourier transform of a function f(x) is defined as 



/(A) := / dxj= 



7(z), 



A e {C : |/(A)| < 00}. 



Assuming /(A) is analytic in an infinite strip in the complex plane A — {A G C : Im(A) G (A_,A+)}, then 
one can recover the function f(x) from /(A) by applying the inverse transform 



f(x) = dA r -J=e^/(A), 



A = A r + iXi, 



where A r , Xi G R. For a review of the g eneralized Fourier transforms a s 



refer the reader to any of the following: 



Bovarchenko and Levendorskiil (2002) 



Xi G (A_,A+), 



they related to Le v y processes, w e 



Lewis! |2001l ): lLiptonl |2002h 



Before presenting explicit solutions uo(t, x) and ui(t,x) we note that 
(A2) 



_1 iXx 

/2¥ L 



b\-k= e iXx 



/2tt 



R 1 „»Aa: 



/27T 



where <j>\ and B\ are given by 

(3.16) cb x = i(j)X - \(o 2 )X 2 + (0 J (e iXz - 1 - iXz) v(dz), 
B x = V 3 (-iX 3 + A 2 ) + U 3 (x 2 J (e z - 1 - z) v(dz) + iX J (e iAz - 1 - iXz^j u(dz) 

(3.17) + V 2 (-A 2 - iX) + U 2 (-iX J (e z - 1 - z)u(dz) + J (e iXz - 1 - iXz^u(dz)J 
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which is valid for any A G C such that <f)\ is well-defined. Those who are familiar with Levy processes will 
recognize <j>\ as the characteristic Levy exponent corresponding to Levy triplet ((7), (c 2 ), (Cjv). 

Theorem 3.1. Assume that <p\, given by (|3.16l) . is analytic in an infinite strip = {A G C : Im(A) G 
(A_ , of the complex plane, which contains the real axis. Assume h(X), the generalized Fourier transform 
of h(x), is analytic in an infinite strip A h = {A G : Im(A) G (A_ , )}. Let A = A r + i\ where A r , A,; G K 
and fix the imaginary component Xi G (A_,AlL). Then the solution uo(t,x) to PIDE Q3.8P with BC (|3.13p is 

u {t,x)= { dA r e*^(A)-i=e a:£! , 

where <f>\ is given by (|3 . 1 6[) . and the solution u\(t,x) to PIDE (|3.11l) with BC (|3 . 14[) is 

u 1 (t,x) = { dX r te^^h(X)B x -k=e lXx . 
Jr V27r 

where B\ is given by (|3.17p . 

Proof. See appendix [XJ □ 



3.3 Accuracy of the pricing approximation 

We have now derived an approximation u £ kj«o+£"i for the price of any European option. However, this 
derivation relied on formal singular perturbation arguments. In what follows, we establish the accuracy of 
our pricing approximation. For our accuracy result, we shall need the following assumption: 

• The payoff function h{x) and all derivatives of h(x) are smooth and bounded. 

Obviously, many common derivatives - e.g., call and put options - do not fit this assumption. To prove 
the accuracy of our pricing approximation for calls and puts would require regularizing the option payoff 



and extending to unbounded functions, as is done in 



Fouaue. Papanicolaou. Sircar, and S0mal (|2003bl) . The 



regularization procedure is beyond the scope of this paper. As such, we limit our analysis to options with 
smooth and bounded payoffs. Our accuracy result is as follows: 

Theorem 3.2. For fixed (t,x,y), there exists a constant C such that for any e < 1 we have 

\u e - (u +em)| < Ce 2 . 

Proof. See appendix [C] □ 



Theorem 13.21 gives us information about how our pricing approximation behaves as e — > 0. In practice e 
is small, but fixed (it does not go to zero). Without knowing what the constant C is in theorem l3.2[ it is 
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difficult to gauge exactly how good the pricing approximation is. As such, in the example provided in section 
21 we will compare the approximate price uq + eu\ of a derivative- asset, calculated using the formulas in 
Theorem 13. 11 to the full price u £ , calculated via Monte Carlo simulation. 



4 Example: Merton jump-diffusion with stochastic volatility and 
jump- intensity 

In this section we provide one specific example within the class of models described in section [2] Specifically, 



we extend the jump-diffusion model of 



Mertonl (|1976T ) to include stochastic volatility and jump-intensity. We 



refer to this class of models as the Extended Merton class or simply ExtMerton. In the Merton jump-diffusion 
model, jumps are log-normally distributed. Thus, we let the measure v be given by 



v(dz) 



1 



: exp 



-{z — m)' 



dz. 



V2^ V 2s 2 
Under this specification, we have 

<7> = -^ 2 }-(0 (> + 4- m -l), 

<t>x = *(7>A - \(a 2 )\ 2 + (C> (>™-^ 2 - iXm - l) , 

B x = V 3 (-iX 3 + A 2 ) + U 3 (A 2 (>+ s2 /2 - l - TO ) + l \ 2a2 /2 _ ! _ iAm )) 

+ V 2 (-A 2 - iX) + U 2 (-iX (e m+s2 ' 2 - 1 - to) + ( e ^m-s 2 A 2 /2 _ 1 _ iAm )) _ 

For a European call option with payoff h(X t ) — (e Xt — e k ) + , the generalized Fourier transform of h(x) is 
given by 

_/c— ikX 



h(X) 



dx —?=e 

V27T 



1 — iXx/x 



Im(A) < -1. 



V2^(iX + X 2 ) 

The values of ((cr 2 ), (C), V3, U3, V2, U2) depend on the particular choice of a(y) and ((y) as well as a specific 
choice for the Y process. In the numerical examples below we let a(y) — —y, P(y) — /3, and A(y) = A so 
that 

dY t = (~Y t - -A A dt+~pdB u 
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and we choose cr(y) = ae v and £(y) = be y . With these choices the invariant distribution of Y under the 

ft 2 

physical measure P is normal Fy ~ N(0, ^-) and we can compute explicitly 

(a 2 ) =aV 2 , 
V 3 = (e? - l) , 

T/ 2 = -/3AaV 2 , 

The implied volatility / corresponding to a European call option with price u is defined implicitly though 

u BS (I) = u, 

where u BS (I) is the price of the call option (with the same strike and maturity) as computed in the Black- 
Scholes framework assuming a volatility of /. In figure Q] we fix the time to maturity at t = 1/10 and we 
plot the Black-Scholes implied volatility smile induced by the approximate price of European calls uq + e u\ 
for different values of e — {0.1, .033, 0.01}. For comparison, we also plot the implied volatility smile induced 
by the full price u £ (computed using Monte Carlo simulation). As expected, as e goes to zero, the implied 
volatility induced by the approximate price uq + e u\ converges to the implied volatility induced by the full 
price u e . 



p 2 

U 2 = -Pkbe^. 



5 Calibration to S&P500 index options 

In this section we calibrate ExtMerton jump-diffusion class discussed in section 2] to the implied volatility 
surface of S&P500 options. For comparison, we also calibra te the classical Merton model and the fast mean- 
reverting stochastic volatility (FMR-SV) class of models of iFouaue. Papanicolaou, and Sircarl (|2000l ) to the 
same set of data. 

In order to formulate the calibration procedure, it will be useful to introduce the following notation 

$ := ({a 2 ),{Q,m,s,V 2 e ,Vi,UI,UI), 

9 := {$ : (a 2 ) > 0, (C) > 0,m G R, a > 0, (V?, Kf , f/f , L/f ) e R 4 }, 



where we have defined Vf := eV{ and Uf :— eUi. Note that the elements of 4> is are the unobservable 
parameters needed to compute the approximate price of an option u + eui in the ExtMerton framework, 
and 8 is the feasible state space of these parameters. Note also that we do not assume a specific value for 
e, a specific volatility process Y, or specific functions: cr(y) or £(?/). In fact, this is one of the main features 
of the class of models considered in this paper. By assuming that the driving factor Y is fast-varying and 
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ergodic, specific choices for (e, Y, o~(y), C(y)) are n °t needed to compute the approximate price uq + eu± 
of an option (or the corresponding implied volatility). For the purposes of calibration and pricing, the 
relevant information about (e, Y, &(y), C(y)) is neatly contained in (c 2 ), (C) and the four group parameters 
{V?,U?,i = l,2}. 

Let I(t, k) be the observed implied volatility of a European call option with time to maturity t and log 
strike k = log-R'. Let I £ (t, k; $) be the implied volatility of a European call option with the same maturity 
and strike as computed in the ExtMerton framework using parameters $ € 6. We formulate the calibration 
problem for the ExtMerton class as a least squares optimization. That is, we seek $* such that 

(5.1) min V (I{t u h) - I 6 (U, h- $)) 2 = V (7(ti, h) - F(t it fcj; #*)) 2 . 

i i 

Here, the sum runs over all pairs (U, ki) for which a call or put is liquidly traded and for which ti > 17 days. 
The calibration procedures for the Merton model and the FMR-SV class are performed in a similar fashion 
by solving (|5.1I) for $ G Q Mer anc | $ g qFMR respectively, where 

e Mer := {$ : (a 2 ) > 0, (C) > 0, m € M, s > 0, (V£, Vi , D|, C/ 3 E ) = 0}, 
6 FMi? := {$ : (a 2 ) > 0, ((C), m, s) = 0, (V 2 E , V?) G K 2 , (f/|, = 0}. 

Note that by requiring (Kf , V^, J7| , f|) = in Mer the effects of stochastic volatility and stochastic jump 
intensity disappear, and the approximate option price in the ExtMerton class uq + e u\ reduces to the Merton 
price uq. Similarly, by requiring that ((C), m, s, C/f , f/f) = in Q FMR ^ the effect of the jumps disappears (the 
effects of stochastic volatility remain), and the approximate option price in the ExtMerton class uo + eu\ 
reduces to the price as computed in the FMR-SV class. 

We perform the calibration procedure for all three frameworks (ExtMerton class, classical Merton model, 
and FMR-SV class) on two dates: October 2, 2006 and May 3, 2010. The dates were chosen to represent 
both the pre-crisis and post-crisis periods. To perform the calibration we use MATLAB's built-in non-linear 
least squares optimizer: lsqnonlin. The fit from the October 2, 2006 calibration (encompassing maturities 
of 47, 75, 166, 257, 355, 446 and 628 days) is plotted in figure [21 The shortest maturity is plotted on the 
top of the figure and the longest maturity at the bottom. For each plot, the units of the horizontal axis are 
log-moneyness to maturity ratio: LMMR := Iog(i4r/S'oe r *). The vertical axis represents implied volatility. 
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The parameters obtained from the October 2, 2006 calibration procedure are listed below: 



Extended Merton : 


(**)- 


= 0.3605 2 , 


(0 = 


= 0.0223, 


m = 


= -0.2226, 


s = 


= 0.4127, 




Vi- 


= -0.0156, 


v 3 s - 


= 0.0001, 


u$ = 


= 0.1304, 


ui-- 


= -0.0481, 


Merton : 


<- 2 > = 


= 0.3605 2 , 


(0- 


= 0.1371, 


m = 


= -0.2750, 


s - 


= 0.1850, 


FMR-SV : 


<- 2 > = 


= 0.3605 2 , 


v 2 -- 


= -0.0767, 


v 3 = 


= -0.0050. 







Extended Merton : 




= 0.1704 2 , 


(0 = 


= 0.0287, 


m = 


= -0.3218, 


s 




Vi- 


= -0.0452, 


Vi- 


= 0.0005, 


ui- 


= 0.0692, 


ui 


Merton : 


(« 2 )- 


= 0.3714 2 , 


(0- 


= 0.0896, 


m = 


= -0.5396, 


s 


FMR-SV : 


<a 2 > = 


= 0.3714 2 , 


v 2 = 


= -0.0547, 


v 3 = 


= -0.0060. 





The fit from the May 3, 2010 calibration (encompassing maturities of 47, 75, 110, 138, 229, 320, 411 and 593 
days) is plotted in figure [31 The parameters obtained from the May 3, 2010 calibration procedure are listed 
below: 

0.7974, 
-0.0528, 
0.3974, 



In both cases, a visual inspection of figures [2] and [3l supports the use of the ExtMerton class over the Merton 
model and the FRM-SV class - especially for high strikes at the shortest maturities. The visual evidence 
is confirmed by the obtained residuals (i.e., the sum in equation (|5.1|) for the optimal <&* obtained in each 
calibration) 

Extended Merton Merton FMR-SV 

October 2, 2006 : 0.0161, 0.0425 0.1306, 

May 3, 2010 : 0.1184, 0.3930 1.0900. 

While, for both dates tested, we find that the ExtMerton class provides a tighter fit to implied volatilities 
than both the Merton model and the FMR-SV class, it is apparent that the calibration for all three models 
performs poorly for the longest maturities. To obtain a tight fit across all maturities, it is likely that multiple 
driving factors - operating on different time scales - are required. We discuss this multiscale extension in 
the following section. 



6 Extension to multiscale stochastic volatility and jump intensity 

The results of this paper can be extended in a straightforward manner to include multiscale stochastic 
volatility and jump intensity. We briefly describe how this may be done. Our intent in this section is not to 
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be rigorous, but rather to give a flavor of the computations involved in this extension. To begin, we modify 
the dynamics of S slightly. Letting S = e x we have 

dXt = 7 (y t , Zt) dt + a(Y u Zt) dW t x + { s dN t (Y t ,Z u ds), X Q = x, 



dZ t = [5 2 c(Z t )-5r(Y u Z t )g(Z t )\dt + 5g(Z t )dW i 





= x, 






Y 


= y, 


> 


(under 


z 


= z. 







Here, Z is a slow-varying factor, in the sense that its infinitesimal generator under P is scaled by S 2 , which 
is assumed to be a small parameter: S 2 << 1. The Brownian motions W x , W y , W z have correlations p xy , 
p xz and p yz (which must be such that the covariance matrix is positive definite), the compensated Poisson 
random measure N(Y, Z, ds) satisfies 

dN t (Y u Z t , ds) = dN t (Y t ,Z t ,ds) - t(Y t , Z t )v(ds)dt, 
E[dNt{Yt,Z t ,ds)\Y u Zt\ = C(Y t ,Z t )u{ds)dt, 

and the drift "f(Y t , Z t ) is given by 

~{{Y u Z t ) = -\o- 2 {Y u Z t ) - aY u Z t ) f (e s - 1 - s)i/(ds). 

Using risk- neutral pricing, the value u £ ' S (t, x, y, z) of a European option in this setting is 

u e < s (t, x, y, z) = E XiViZ [h(X t )\ , h(x) := H(e x ). 

From the Kolmogorov backward equation, the function u e ' S satisfies the following PIDE and BC 

(6.1) {-d t +A E ' S ) u £ ^ & = 0, u E ' S (0,x,y,z) = h{x), 

where the partial integro-differential operator A £ ' 5 is the generator of {X, Y, Z), The operator A e ' 5 has the 



following form 



A e ' & = — A + -Ai + A 2 + -M 3 + SM 1 + S 2 M 2 . 



Terms containing 5 in (|6.ip are small in the small-<5 limit, giving rise to a regular perturbation. Thus, (|6.1[) 
has the form of a combined singular-regular pe rturbation about the 0(1) operator (—dt +A%). Following 



Fouque. Papanicolaou. Sircar, and Solna 



(|201lh we seek a solution u e ' of the form 



oo oo 
7i—0 m— 
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Our goal is to find an approximation u £ ' S — uq^q + £ Ui,o + & u o,i + 0(e 2 + S 2 ). A formal asymptotic analysis 
yields the following PIDEs for uq.o, ui.o an d 1*0,1 

0(1): {-d t + (A 2 ))u ,o = 0, uo tQ (0,x,z) = h{x), 

0(e): (-dt + (A 2 ))ui, = -Su 0j o, «i,o(0, s, 2) = 0, 

0(<5): (-9t + <X 2 »uo,i = -<M 1 >u ,o, «o,i(0, s, z) = 0, 

where, as in section |3"7TI the y-dependence has disappeared from uo,o ; ui,o and Uo,i- The operators (.A2), 2? 
and (Mi) are given by 



(A 2 ) = (-,(■, z))d x + ha 2 (-,z))dt x + (((;z)) / (e'«"-l-*a s ) !/(*»), 
z Jr 

® = ^ 3 (z) - d 2 xx ) + U 3 (z) (-j^(e s -l- s)v{ds)dl x +J(o a -l- sd x )d x v(ds) 

+ V 2 (z) (d 2 xx - d x ) + U 2 (z) f-J^(e s -l- S y(ds)d x + J (9 S - 1 - sd B y(ds)j 
(Mi) - -g(z)(r(;z))d z +g(z)p xz (a(- lZ ))dl z , 
where the z-dependent parameters {V 3 {z), U 3 {z), V 2 (z), U 2 (z)) are 

V 3 (z) = -^L((3(.)a(;z)d y r ] (;z)), U 3 (z) = -p xy {f3(-)a(;z)d y Z(;z)), 



V 2 (z) = -(j3(-)A(;z)d yV (;z)), U 2 (z) = (/3(-)A(;z)d v £(;z)). 



The expressions for u q and ui are analogous to those given for uq and U\ in Theorem 13. II An expression 
for uo ; i is obtained using the theory generalized Fourier transforms 

v(s,X,z) := J dx^=e~ lXx (Mi)u Q ^(s,x,z) 

u 0t i(t,x,z)= f d\ rT ±=e lXx f dse (t - s) ^ {z) v{s,\z). 







where, as in Theorem 13.11 the Fourier variable satisfies A = A r + i\ where Ai 6 (A^ , Note, care must 
be taken when computing (Mi)uo j o as both terms in (Mi) contain the operator d z and wo,o depends on z 
through both (tr 2 (-,z)) and {(,(■, z)). A careful computation shows that u .i is linear in the following four 
parameters 

V 1 (z)=g(z)p xz (a(- 7 z))d z (a 2 (-,z)) 7 V (z) = -g(z)(T(; z))d z {o 2 {; z)), 

U x {z) = g(z) Pxz (a(;z))d z (a;z)), U (z) = -g(z)(T(-, z))d z {((; *))■ 
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Finally, the accuracy of the multiscale pricing approximation uq,q + eiti.o + Sug,i is as follows: for fixed 
(t, x, y, z) there exists a constant C such that for any e < 1, 8 < 1 we have 



u e ' s - (u 0;0 + £u 1>0 + 5u ,i)\ < C{e z + 8 Z ). 



The p roof of this error bound is analogous to the proof found in chapter 4 of 
POIII) . 



Fouaue, Papanicolaou. Sircar, and Solna 



7 Conclusion 

In this paper, we have introduced a class of exponential Levy-type models in which the volatility and jump- 
intensity are driven stochastically by two factors - one fast-varying and one slow-varying. Using techniques 
from the theory of generalized Fourier transforms, singular and regular perturbation theory we have derived 
a general formula for the approximate price of a European-style derivative. Furthermore, we have quantified 
the accuracy of our pricing approximation both theoretically (see Theorem 13. 2[) and numerically (see figure 
[I}. We hope this work motivates further research into exponential Levy- type models. A possible extension of 
this paper would be to allow the jump distribution (rather than just the jump intensity) to vary stochastically 
in time. 
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A Proof the Theorem 13.11 

We wish to solve PIDE (13.81) with BC (|3.13[) . Throughout this appendix, we always assume A = A,. + iXi 
where Xi G (A_, A'j:). Under the assumptions of Theorem 13. II we have 

d t u (t,x) = (A 2 )u (t,x) =■ d t [dx^=e- iXx u {t,x) = jf dx -j= e - iXx (A 2 )u {t, x) 

dxu (t,x)(A 2 )*^e- lXx 

= 0A [ dxu (t,x)^=e- lXx 

(A.l) =► ftuo(t,A) = ^A«o(t,A), 

(A.2) u (Q,x) = ft(x) u (0, A) = ft (A). 

Note that (A2}*, the formal adjoint of (.A2}, is given by making the replacements d x — > —d x and Z — > 9- z 
in (|3.15|l . Using (|A.1|) and (|A.2|) one deduces 

u Q (t,X) =e t,px h(X) => u (t,x) = { dX r e tc/> ^h(X)^=e iXx . 

To solve PIDE (l3TTTj) with BC (|3"H)) we first note that 

■Bu Q (t,x) = J d/i r e t +»h(ji)$^-e i '" = J d^e^h^B^e 1 ^, 

where /1 = \i r + iXi. Thus, 

/ dx^e- iX *Vu (t,x)= f dHre»»h(M)Bp f dx j=e~ iX * j=e^ 

G?/i r e t ^ M /i(/j,)B A1 (5(/i r — A r ) 
= e^ft(A)S A , 

where we have used ^ / R dx e ~KK-Hr)x _ s ^ r - Mr ). Therefore multiplying PIDE (|3~TT|) and BC (j3~T4| by 
_|^ e -iAa; an( j integrating with respect to x one finds 

(-9* + 0a)wi(*,A) = -e t ^h(X)B x and u x (0,A)=0. 

Solving for U\{t, A) we have 

Ui(t,A) =te t4 ' x 7i(X)Bx => Ul (t,x) = [ dX r te t4 "<h(X)Bx-k=e iXx , 

which completes the proof. 
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B Existence of classical solutions Cauchy problems 



A classical existence and uniqueness result for second o rder integro-differential 



coefficients is established in Theorem 3.2, Chapter 3 of 



Cauch y problems with regular 



Bensoussan and Lionsl (|1984l ). Below, we present the 



results of this Theorem in a form which is convenient for our purposes. We introduce the following function 
spaces: 

• W 2 - p {M. n ) is the space of functions u such that u <E L p (M. n ), d l u € L p (R n ) for all i and 9?it E L P (R™) 
for all 



• LP(0, T; W 2 ' P (W 1 )) is the space of functions u such that 

• L p (0,T; L P (R™)) is the space of functions u such that 



\LP( 



LP((0,T)) 
< CO. 



< 00. 



Lp((0,T)) 

W 2 ' 1,P ((Q,T) x R") is the space of functions u G L p (0, T; W 2 ' p (M. n )) such that d t u € L p (0, T; L P (W 1 )). 



The above norms are the usual L p (0) and W k ' p (fi) norms (below, a is a multi- index) 



IMIlp(Q) = ( I q \u(x)\ p dx 



i/p 



i/p 



,|a|<2 



Theorem B.l. Let A be given by 

n n „ 

A =Y1 ai i ( x ) di i + H b i(x)di + c{x) / v{dz) (0 Z - 1 - z ■ V) , 

i,j=l i=l 

where x € R". Assume 

aij,bi,c € L°°(R") and V" CLijXiXj > a\x\ 2 , Vrr € R", a > 0. 

ITien /or 3 € L 2 (0, T; L P (R™)) and ft G M /2,P (R"), 2 < p < oo, i/ie problem 

(-d t +A)u = g, u{0,x) = h(x), u e W 2 ' 1,p ((0,T) x R"), 

/ias a unique classical solution. Furthermore, if g,h are bounded, then u is bounded. 



C Proof of accuracy 

Before establishing our main accuracy result - Theorem 13.21 - we shall need the following lemma. 

Lemma C.l. Suppose J{y) is at most poly nomially growing. Then, under the assumptions of section \2.1\ 

for every y and s < t, there exists a positive constant C < oo such that for any e < 1, we have the following 
inequality 

E y [\J(Y s )\]<C. 
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Proof of Lemma [C.l[ See Lemma 4.9 of lFouaue. Papanicolaou. Sircar, and Solnal ([20111 ). 

We are now in a position to prove Theorem 13.21 We begin by defining a remainder term R £ by 



□ 



u £ = u + e u\ + e u 2 + £ u 3 + R £ 



From the equations satisfied by uo, ui, u 2 and u 3 in section f3.ll one derives 



(C.l) 



(-d t +A £ )R £ = -e 2 F £ , 



where we have denned 



F £ (t, x, y) := (Aiu 3 + A 2 u 2 ) + eA 2 u 3 , 



R £ (0,x,y) = -e 2 G £ (x,y). 



G £ {x, y) := u 2 (0, x, y) + eu 3 (0, x, y). 



Now, fix some finite T. For t € [0, T] we note that F £ and G £ are smooth functions of t,x,y that are, 
for £ < 1, bounded by smooth functions of t,x,y independent of e, uniformly bounded in t,x and at most 
polynomially growing in y (see lFouaue. Papanicolaou. Sircar, and Solnal ([20111 ). p. 143). This assertion can 
be deduced as follows. First, since h is smooth and bounded (by assumption), then by Theorem IB. 11 the 
function uq is bounded. Likewise, since 9™ commutes with (A 2 ), and since all derivatives <9™/i are bounded 
(by assumption), then by Theorem lB.il one deduces that 8"uq is bounded for every n. Now, we note that 
u 2 and U3 are completely characterized by uq and ui through Poisson equations (13.51) and p.6|) . From (|3.10[) 
one finds that u 2 is given by 



"2 



- Q»? (9 2 xx - d x ) - ^ (e z - 1 - z) v(dz)d x +ijj0z 



— 1 — zd x ) v{dz) I u + D 



where D(t,x) is an arbitrary function, which we choose to be zero (the choice of D does not effect our 
accuracy result). Because all derivatives d x uo are bounded, one deduces that u 2 is uniformly bounded in 
t, x. The y-dependence in u 2 enters only through the solutions r\ and £ of Poisson equations (|2 . 2[) . which we 
assumed are at most polynomially growing in y. Likewise, derivatives d x u 2 and <9™ u 2 are uniformly bounded 
in t 1 x and at most polynomially growing in y. A similar argument shows that u 3 and its derivatives d"u 3 
and dyU 3 are uniformly bounded in t, x and at most polynomially growing in y. Now, note that F £ also 
has y-dependence through the operators A\ and A 2 . These operators depend on y through the coefficients 
a, A, £, which are all bounded by assumption, and through /3, which is at most polynomially growing in y 
(again, by assumption). Thus, as claimed, F £ and G £ (as well as u 2 and u 3 ) are uniformly bounded in t,x 
and at most polynomially growing in y. And, for e < 1, the functions F £ and G £ are bounded by functions 
of t, x, y, which, independent of £, are uniformly bounded in t, x and at most polynomially growing in y. 
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Now, note that the solution R £ of (|C.1[) has the following Feynman-Kac representation 



(C.2) 



R £ (t,x,y) = e 2 E : 



■x,y 



G £ (X u Y t )+ [ F £ (s,X s ,Y s )ds 



The validity of the Feynman-Kac representation is guaranteed because (i) the expectation (|C.2[) is clearly 
finite (due to the uniform bound in s < t and x for F £ and in x for G £ , the polynomial growth bound 
in y for F £ and G £ , and the existence of moments E[|y| fc ] of any order), and (ii) the function R £ := 
u £ — (uq + eu\ + e 2 it2 + £ 3 u 3 ) is smooth by construction. 

Using representation (|C.2[) for R £ and Lemma ICTTl there exists a constant C\ < oo such that 



for some constant C < oo. This concludes the proof of accuracy. 
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Figure 1: Using the model described in sectional we plot the implied volatility induced by the price of 
European call option as a function of the strike price K. In each plot, the dashed blue line corresponds 
to the implied volatility induced by the full price u £ (computed via Monte Carlo simulation) and the solid 
green line corresponds to the implied volatility induced by our approximation uq + eu\. For all plots we use 
the following parameter values: t = 1/10, e x = 50, m = —0.2, s = 0.2, p = —0.7, a = 0.2, b — 1.5, (3 = 1.0 
and A = 0.25. 
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Figure 2: Implied volatility fit to S&P500 options from October 2, 2006. 
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Figure 3: Implied volatility fit to S&P500 options from May 3, 2010. 
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